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ABSTRACT 

We present a new catalogue of galaxy triplets derived from the Sloan Digital Sky Survey 
Data Release 7. The identification of systems was performed considering galaxies brighter 
than Mr = —20.5 and imposing constraints over the projected distances, radial velocity dif- 
ferences of neighbouring galaxies and isolation. To improve the identification of triplets we 
employed a data pixelization scheme, which allows to handle large amounts of data as in 
the SDSS photometric survey. Using spectroscopic and photometric data in the redshift range 
0.01 ^ 2 ^ 0.40 we obtain 5901 triplet candidates. We have used a mock catalogue to analyse 
the completeness and contamination of our methods. The results show a high level of com- 
pleteness 80%) and low contamination 5%). By using photometric and spectroscopic 
data we have also addressed the effects of fiber collisions in the spectroscopic sample. We 
have defined an isolation criterion considering the distance of the triplet brightest galaxy to 
closest neighbour cluster, to describe a global environment, as well as the galaxies within a 
fixed aperture, around the triplet brightest galaxy, to measure the local environment. The final 
catalogue comprises 1092 isolated triplets of galaxies in the redshift range 0.01 ^ z ^ 0.40. 
Our results show that photometric redshifts provide very useful information, allowing to com- 
plete the sample of nearby systems whose detection is affected by fiber collisions, as well 
as extending the detection of triplets to large distances, where spectroscopic redshifts are not 
available. 
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1 INTRODUCTION 

In a hierarchical scenario, the interactions affects general proper- 
ties of the galaxies such as stor formation rate, nuclear activity and 

morphology dAlonso et al.l2004l.l2006h 

Al though galaxy pairs jLambas et alj 1200 3l : lAlonso et al.l 

|2004 
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galaxy pairs 
and gr oups 



with four or more members 

( iMerchan & Zandivarea |2005|) have been the subject of sev- 
eral papers in the literature, there are far few works on triplets 
of galaxies, a link between pairs and groups and a potentially 
interesting laboratory to investigate several process associated to 
galaxy evolution. 

Traditionally, galaxy groups selection have excluded triplets 
of galaxies because these systems were not considered as bonafide 
physical structures, i.e., dynamically isolated systems. Neverthe- 



E-mail: aomill@oac.uncor.edu 



less. lHemandez-Toledo et"aL have recently analysed a sam- 

pl e of galaxy triplets der ived from the catalogue of isolated triplets 
of lKarachentseval jl973h finding signatures of physical interactions 
(tidal bridges, excess of barred galaxies), that implies that most of 
these objects are indeed real physical structures. 

Pioneering works on t he identification of triple ts of 
galaxies were performed by iKarachentseva et si\ ( 1 1979 ^ and 
[Karachentsev et al. ( 1988). These authors present a list of 84 North- 
ern isolated galaxy triplets with members with component ap- 
parent magnitudes brighter than 15.7, selected by visual inspec- 
tion of Palomar Sky Survey prints. They found that about 24% 
of the members are elliptical and lenticular galaxies, while 76% 
are spirals and irregulars. Numer ous studies have been made in 
the fo llow up of this catalogue. iKarachentsev & Karachentseval 
( Il98lh measured radial velocities for the isolated triplets and found 
that 5% of the systems are spurious (AVij > 500km s~^), 
31% consist of a pair of galaxies with nearly equal radial ve- 
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locities (AVij < 500km s^^) and a projected third component 
and 64% are considered physical triplets (AVij < 500km s~^). 
Velocity dispersion, diameters, integrated luminosity, virial mass, 
and the virial mas s to luminosity ratio were estima t ed fo r 
the physical triplets jKarachentseva & Karachentsevll 19821 1 1983h . 
The spatial configuration, dynamics and dark matter c ontent of 
these systems was studied by several authors ( Karache ntsev et al. 
1 990 ; Chernin & Mi kkolj|l99ll : lAnosova et a"illl992l ; IZheng et al. 
1993; Aceves 2000|). For the Sou thern sky ((5 < -3°) 
Karachentseva & Karachent"sevl j2000t) selected 76 isolated triple 
systems of galaxies using the ESO/SERC and POSS-I sky surveys, 
with the same criteria defined for t the Northern systems. 

iTrofimov & CheminI j 19951) compiled a list of 108 triplets 
of galaxies sel e cted from two group cat alogues, those of 
iGeller & Huchral jl983h (Northern sky) and iMaia et all h989l) 
(Southern sky). 38 of these systems are considered isolated. The au- 
thors call these systems "wide triplets" in contrast to the compact 
triplets (with smaller mean projected harmonic separation) stud- 
ied by iKarachentseva et al.l The main assumed difference 
between wide and compact systems, regardless their difference in 
size, is that galaxies in compact triplets are supposed to have in- 
teracted many times and are probably in equilibrium within the 
system. In contrast, wide triplets would be dynamically young and 
probably far from virialization. 

A high-or d er 3 D Voronoi tessellation method was employed 
by lElviv et aU (2009) for identifying single galaxies, pairs and 
triplets on a volume-limited sample of galaxies from the Sloan Dig- 
ital Sky Survey Data Release 5. Nevertheless, a catalogue of triplet 
of galaxies has not been compiled for more recent SDSS releases. 

In this work we present a catalogue of triplets of bright galax- 
ies {Mr < —20.5) selected from a volume limited sample extracted 
from SDSS-DR7 in the redshift range 0.01 2 0.4. We used 
data with both sp ectroscopic redshift measurements and photomet- 
ric redshifts from lO'Mill et al.l ( I2OIII) . 

This paper is organized as follows. In Section 2 we describe 
the galaxy samples used in this work. Section 3 describes the algo- 
rithm developed for the detection of triplet of galaxies candidates. 
The implementation of this algorithm at low redshifts is described 
in Section 4. In this section we also perform a completeness and 
contamination test to our algorithm and analyse the incompleteness 
effect due to fiber collisions. Section 5 describes the extension of 
our results to higher redshifts. In Section 6 we describe the isolation 
criteria employed to define the final sample of triplet of galaxies. 
Finally, Section 7 summarized the results obtained in this work. 

Throughout this paper we adopt a cosmological model char- 
acterised by the parameters Q,rn = 0.25, SIa ~ 0.75 and Hq = 
70km s"^ Mpc"^ 



2 SPECTROSCOPIC AND PHOTOMETRIC GALAXY 
SAMPLES 

The samples of galaxies used in this work were drawn from 
the Data Release 7 of Sloa n Digital Sky Su rvey (SDSS-DR7, 
lAbazaiian et"ai]|2009l) . SDSS jYork et alj[2000b has mapped more 
than one-quarter of the entire sky, performing photometry and spec- 
troscopy for galaxies, quasars and stars. SDSS-DR7 is the seventh 
major data release, corresponding to the completion of the sur- 
vey SDSS-II. It comprises 11.663 sq. deg. of imaging data, with 
an increment of ~ 2000 sq. deg., over the previous data release, 
mostly in regions of low Galactic latitude. SDSS-DR7 provides 
imaging data for 357 million distinct objects in five bands, iigriz. 



as well as spectroscopy over ~ tt steradians in the North Galactic 
cap and 250 square degrees in the South Galactic cap. The aver- 
age wavelengths corresponding to the five broad bands ar e 3551, 
4686, 6165, 7481, and 8931 A (Fukugita et al. 1996; Ho gg et al. I 
I2OOI; Smith et al. 20(^. For details regarding the SD SS camera 
see Gunn e t al., t I998i) ; for astrometric calibrations see IPier et all 
12003). The survey has spectroscopy over 9380 sq. deg.; the spec- 
troscopy is now complete over a large contiguous area of the North- 
ern Galactic Cap, closing the gap that was present in previous data 
releases. 

In this work we employed spectroscopic and photometric 
data extracted from SDSS-DR7. The spectr oscop ic data were de- 
rived from the Main Galaxy Sample (MGS; IStrauss et al. I ( i2002h ) 
obtained from the fits files at the SDSS home pagqj. k- 
corrections for this sa mple were calcula t ed usi ng the software 
k-correct_v4 . 2 of iBlanton & Roweij ilOOH ). The photomet- 
ric data_were_derived from the photometric catalogue constructed 
bv lO'Mill et al.l Jioi llPl These authors compute photometric red- 
shift and k-correction for the photometric data of the SDSS-DR7. 
The rms of the photometric redshift is (Jphot ~ 0.0227 and k- 
corrections were obtained through joint parametrisation of redshift 
and reference frame (at 2 — 0.1) (g — r) colour. For both data 
sets, k-corrected absolute magnitudes were calculated from Pet- 
rosian apparent magnitudes converted to the AB system. 

In order to use these spectroscopic and photometric data, it is 
necessary to analyse the completeness and reliability of these sam- 
ples. To explore this issue we have calculated the spatial number 
density of galaxies brighter than Mr = —20.5 for both spectro- 
scopic and photometric data, taking into account the apparent mag- 
nitude limit r = 17.77 and r = 21.5, respectively. In Figure[T]it 
can be appreciated that the spatial number density of galaxy for, 
both, spectroscopic and photometric samples, shows a flat trend up 
to z ^ 0.14 and z ~ 0.4 respectively. The spatial number density 
of galaxies for the photometric data shown in Figure 1 corresponds 
to a random sample with ~ 0.001% of the complete photometric 
sample. 

Taking into account these considerations we derived four 
galaxy data sets from SDSS-DR7: 

The first sample, hereafter SI, comprises galaxies brighter 
than r — 17.77 with spectroscopic redshift in the range 0.01 ^ 
-Zspec < 0.14. This sample was employed for the detection of triple 
systems of galaxies candidates at low redshift. 

In order to analyse the effect of photometric errors over the 
radial velocity differences AV for the photometric data, we have 
selected a sample that comprises galaxies brighter than r — 17.77 
with spectroscopic measurements in the redshift range 0.07 ^ 
Zspcc ^ 0.12 for which we have obtained photometric redshifts. 
We will refer to this sample hereafter as the S2 sample. 

We have also extracted a sample of galaxies brighter than r — 
17.77, that comprises galaxies with spectroscopic measurements in 
the redshift range 0.09 ^ ^spoc ^ 0.11 and galaxies with photo- 
metric redshifts in the range 0.09 — Uphot ^ ^phot ^ 0.11-1- (Jphot 
(hereafter S3). This sample was employed in the analysis of com- 
pleteness due to fiber collision, described in Section 4. 

Finally, for the detection of triplets of galaxies candidates at 
intermediate redshift we employed a sam ple derived from the pho- 
tometric catalogue o fl O'MilletallJioT l|). We have selected a sam- 
ple that comprises galaxies brighter than r — 21.5, in the redshift 
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Figure 1. Spatial number density of galaxies brighter than Mr = —20.5. Left: Spectroscopic data with apparent magnitude r < 17.77. Right: Random 
sample of ~ 0.001% of the photometric sample with appai'ent magnitude r < 21.5. En'or bars represent Icr uncertainties calculated by bootstrap. The 
vertical line shows the limits chosen for the detection of the triple candidate systems. 



range 0.09 ^ Zphot ^ 0.5. We will refer to this sample as S4 
hereafter. 

Table 1 summarises the main characteristics of the different 
samples used in this work. 



3 BUILDING THE CATALOGUE 

In this section we present the algorithm adopted to identify galaxy 
triplet candidates from the spectroscopic and photometric volume 
limited samples of SDSS-DR7. The algorithm searches for triplets 
by selecting galaxies within a given projected distance (rp) and 
radi al velocity diff erence (AV) of another galaxy. 

iLambas etai] (|2003) identified pairs of galaxies using = 
100 kpc and AV = 350km s"\ iMerchan & ZandivarezI 
l l2005h identified galaxy groups in the SDSS third data release 
with a mean redshift of 0.1 and a median velocity dispersion of 
230km s~^. These authors employed a fiends-of-friends algorithm 
with a transverse linking length corresponding to an overdensity of 
80 and a line-of-sight linking length of AV = 200km s~^. For 
the detection of triple galaxy systems, we have considered these 
restrictions in the selection of Vp and AV constrains. 



3.1 Data Pixelization 

For the purpose of a suitable selection of the systems, as well as 
to efficiently reduce the computing time, we have used a pixeliza- 
tion of SDSS data. Indeed, data pixelization is an important tool in 
managing large amounts of information such as that contained in 
the SDSS photometric catalogue. 

In order to perform a pixelization of SDSS data, we have em- 
ployed routines of the Hierar chical Equal Area Pixelization iso- 
Latitude software (HEALPIX, fG6rski et allboOSl) . HEALPIX was 
developed for data processing and analysis of the observations of 



the cosmic microwave background (CMB). This pixelization pro- 
cedure consists in an equal area isolatitud partition of the sphere, 
which results into a versatile structure for data analysis. The reso- 
lution base comprises 12 pixels formed by three rings around the 
poles and the equator. For higher resolutions, each pixel is subdi- 
vided into four equal-area pixels of smaller size. The grid resolu- 
tion is set by the N^ide parameter, which defines the number of 
divisions along the side of a pixel, required to achieve the high- 
resolution partition desired. All pixel centres are equidistant in az- 
imuth for each ring and are placed on rings of constant latitude. 
All rings located in the equatorial zone are divided in the same 
number of pixels, A^e, = iNside- The rings located in the polar 
cap region contain a varying number of pixels that increase by one 
pixel within each quadrant, with increasing distance from the poles. 
The resulting map contains Npix = VIN'^^^^ pixels of equal area 

HEALPIX has a library of computational algorithms and vi- 
sualisation software that allows fast scientific applications on dis- 
cretized maps created from large amounts of data as is the SDSS 
survey. 

We have employed HEALPIX routines to find the index of 
all pixels within a radial angular distance from a defined center. 
The routines used were: ang2pix, which generates a linking list 
for galaxies in each pixel; ang2vec, that converts angular to Carte- 
sian coordinates, and, finally, query jdisc, that identifies the central 
galaxy pixel and the adjacent pixels inside the search radius, pre- 
viously defined. The resolution (Nside) was set to 512 at redshift 
smaller than 0.14 (SI, S2 and S3 samples) and 2048 at redshift 
greater than 0.14 (S4 sample). The choice of different values of 
Nside is due to the density difference between the lower and inter- 
mediate redshift samples. 

For each selected center we define a search radius that covers 
approximately 600 kpc at the smallest redshift of each sam- 
ple, ensuring the search in every pixel adjacent to the center. The 
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Table 1. Description of the samples used in this work. 



sample name 


Description 


number of objects 


SI 


MGS spectroscopic data. 

Used for the identification of triplets of galaxies candidates at low redshift (0.01 ^ 2:spcc < 0.14) 


684939 


S2 


MGS with photometric redshift (0.07 ^ ^spec s£ 0.12). 
Used in section 4.3 


248210 


S3 


MGS and photometric data. 
Used in the fiber collision analysis 


80851 and 860021 


S4 


Photometric catalogue with redshifts and k-correction from O'Mill et al. (2011) 

Used for the identification of triplets of galaxies candidates at intermediate redshift (0.09 ^ Zphot 


66390174 

^ 0.5) 
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Figure 2. Distribution of AVrma for triplet candidates identified in the spectroscopic sample SI. Red (continuous) line: Tp < 100 kpc and AV < 
500km s~^. Blue (dashed) line: systems with Tp < 200 kpc and AV < 700km s~l. 



adopted search radii were: 0.88 deg for S 1 , S2 and S3 samples and 
0.068 deg for the S4 sample. 



3.2 The algorithm 

In this section we describe the successive steps of the algorithm 
developed for the identification of triplet of galaxies candidates. 

Step 1: The algorithm begins with a selection of central 
galaxy candidates. The pixel of each candidate is identified. Then 
a search for galaxies belonging to adjacent pixels, enclosed within 
the search radius, is performed. In the first run, any galaxy is con- 
sidered a central galaxy candidate. 

Step 2: To identify triplet systems, we have considered as 
neighbours all those galaxies that satisfy the constraints in pro- 
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jected distance and radial velocity difference, ie. Vp < Tpmax and 

AV < AVrnax- 

Step 3: We consider a triplet candidate when there are only 
three galaxies within the constraints in Step 2. 

Step 4: Then, we assign the brightest galaxy to the center of 
each system. This generate two types of multiple identifications; 

• The center is assigned to more than one system indepen- 
dently identified. This is because the brightest galaxy is a neighbour 
galaxy in more than one system. 

• A system will be identified more than once if different cen- 
tres satisfy the neighbourhood condition, ie have Vp < rpmax and 

AV < AVmax- 

In order to avoid double identifications, repeated centres were re- 
moved from the sample and close galaxy pair centres were consid- 
ered as a single system by taking the brightest galaxy as a central 
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galaxy. This procedure generates a new sample of central galaxies. 
With these centres as input sample, the algorithm returns to Step 1 . 

Step 5: It may happen that the new central galaxies may be 
located at a distance greater than Vpmax or have a radial velocity 
difference greater than AVmax with respect to one of the selected 
neighbours. To retrieve these systems, we consider a second set of 
constraints, relaxing restrictions over velocity difference and pro- 
jected distance, and then returning to Step 2. 



4 IMPLEMENTATION OF THE ALGORITHM AT LOW 
REDSHIFT (0.01 2spoc < 0.14) 

The aim of this work is the identification of triplets of luminous 
galaxies, defined as systems close in projected distance and radial 
velocities. Building of a catalogue of triplets with spectroscopic 
measurements is less subject to contamination than using photo- 
metric redshifts given the small uncertainties. So, in principle, us- 
ing spectroscopic data results into a more reliable identification of 
these objects. 

4.1 Spectroscopic data 

In order to identify triplets of galaxies candidates in the S 1 sample, 
we have run the algorithm considering the following restrictions 
for the selection of neighbouring galaxies: Vp < 100 kpc and 
AV < 500km s'^ 

Following the steps described in the previous section, we have 
assigned the center to the brightest galaxy of the system. After 
this recentering, we have considered a second group of constrains 
rp < 200 kpc and AV < 700km s~^, in order to retrieve the 
systems lost by the procedure. 

Therefore, a total of 273 triple system candidates were ob- 
tained from the SI sample, 164 of them with the constraints 
rp < 100 h-^ kpc and AV < 500km s"^ and 109 with Vp < 
200 /i"^ kpc and AV < 700km s'V 

Figure |2] shows the distribution of AVrms, where AVrms 
is the root mean square radial velocity difference for the sys- 
tems. Red lines correspond to the triplet candidates obtained whit 
rp < 100 kpc and AV < 500km and in blue lines 
we present the systems obtained with rp < 200 h^^ kpc and 
AV < 700km s^^. Both distributions present a similar behaviour. 

4.2 Completeness and contamination tests 

In order to test the algorithms, we performed a completeness and 
contamination test employing a SDSS mock catalogue derived 
from the semi-analytic models of Croton et al. (2006), obtained 
from the Millennium Simulation ( Springel e t alj 120051) outputs. 
The spatial resolution of this simulation is suitable for the imple- 
mentation of our algorithm. The luminosity function of the semi- 
analytic model is consistent with observations of galaxies at low 
redshift. The mock catalogue comprise 514914 galaxies brighter 
than r — 17.77, where r correspond to the r — band magnitude 
of the SDSS. This mock catalogue contains: right ascension, dec- 
lination, apparent magnitude in r ~ band, halo mass, redshift and 
peculiar velocity. 

For the completeness and contamination test we have em- 
ployed a volume limited mock catalogue considering galaxies 
brighter than Mr = -20.5 in the redshift range 0.01 ^ z < 0.14. 

The completeness is defined as Nrec/ Ntot, where Nrec is the 
number of triplet candidates identified with our algorithm that are 



matched to the systems found in real space (recovered triplets), and 
Ntot is the number of systems selected with real distances (real 
triplets) within the same halo. We find that the completeness 
achieved by our algorithm is ~ 80%. 

The contamination is defined as Nsp/Ndet, where Ndet is the 
number of triplet candidates identified with rp and AV^. Nsp is 
number of systems detected whose member galaxies do not be- 
long to the same halo (spurious triplets). The contamination for our 
triplet sample is very small, ~ 5%. 

Figure[3]shows the redshift distribution for the real, recovered 
and spurious triplet of galaxies, as well as the dependence of the 
completeness and contamination rates with redshift. 

4.3 Analysis of photometric data in the spectroscopic 
redshift range 

The major source of contamination in the detection of the photo- 
metric systems is due to the error in Zphot, since the uncertainties 
in photometric redshifts affects the choice of AV. To evaluate this 
point, we employed the S2 sample. We have run the algorithm for 
different values of AV, corresponding to laphot * c, 1.5aphot * c 
and 2aphot * c, where Uphot is the mean photometric redshift error 
{lyphot ^ 0.0227, as previously mentioned in Section 2) and c is 
the speed of light. We have then examined the distribution of the 
values of the rms spectroscopic velocity differences ctav, consid- 
ering the spectroscopic redshifts of the triplet candidate members 
identified photometrically. 

For the S2 sample we identified 80 triplet candidates that 
satisfy the constraints described en section 4.1. Using photomet- 
ric data we recover 71 of these systems considering laphot and 
l-5o"phot for AV. For 2Gphot we recovered all the spectroscopic 
triplets. 

In figure|4]we show the gav distribution of triplet candidates 
identified in S2 sample. In different lines and colours we show 
the systems identified with laphot in red, l.^aphot in magenta, 
and 2aphot in blue line. From the figure it can be appreciated that 
we have a large level of contamination when using l.^cTphot and 
2aphot- Therefore, we conclude that the laphot interval provides a 
suitable compromise between high completeness and low contam- 
ination. 

4.4 Incompleteness due to fiber collisions 

The SDSS spectrograph uses fibers manually connected to plates 
in the telescope's focal plane. T hese fibers are mapped through a 
mosaic algorithm jSlanton et al.f 2003) that optimises the observa- 
tion of lar ge-scale structures. Two fibers can not be placed closer 
than 55" ( iStrauss et al. Il20()2h . so for two objects with the same 
priority (such as two MGS galaxies) and whose centres are closer 
than 55", the algorithm selects at random the galaxy which will 
be observed spectroscopically. There are regions where the plates 
overlap (about 30% of the mosaic regions), in which both objects 
may be observed. 

Due to fiber collision the spectroscopic sample is affected by 
incompleteness. The magnitude limit of spectroscopic objects is 
r = 17.77, but not all galaxies brighter than this limit were ob- 
served. This issue becomes more important when analysing com- 
pact objects. 

Tp = r X cos{9) X do X d'^, where r, and <I> are the spherical coor- 
dinates with origin in the triplets center. 
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Figure 3. Completeness and contamination analysis. Left panel: Real, recovered and spurious systems redshift distributions. Right panel: Completeness (solid 
line) and contamination (dashed line) as function of redshift. 
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Figure 4. ct/wa spectroscopic distribution for photometric triplet candidates using S2 sample. Black (short dashed line): spectroscopic systems. Red (dotted), 
magenta (dot-short dashed line) and blue (short dashed-long dashed line): photometric systems obtained considering Icp^g^ * c, LScpjiot * c and 2(Tpf^^^ * c 
velocity intervals, respectively. 
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Figu re 5. Spatial numbe r density of galaxies as a function redshift for the MGS (red dots) and for a random sample with 0.001% of the photometric sample 
from lO'Mill et alj j201 ih (blue dots) both with r < 17.77. The error bars represent la uncertainties calculated using bootstrap. For both samples the galaxy 
density for Zspec < 0.14 is nearly flat. The small and large regions correspond to redshift intervals employed in the analysis of fiber collision (S3 sample). 



IZehavi et alj ( I2OO2I) found that approximately 6% of the tar- 
get objects are no t assigned to a fiber due to the 55" restric- 
tion. lAlonso et alj l|2006) have identified pairs of galaxies in the 
MGS of SDSS second data release and analysed galaxy pairs loss 
due to fiber collision. These authors found that approximately 
5% of the pairs were not recovered due to this effect. Recently 
ICosta-Duarte et al. I jlOllh found that ~ 40% of galaxies with 
r < 17.77 have spectroscopic measurements. 

In figure |5] we show the spatial galaxy density as a function 
of redshift for the MGS (red dots) and for a random sampl e with 
0.001% of the photometric data from lO'Mill et al.l ( 1201 ll) (blue 
dots) both with r < 17.77. It can be appreciated a nearly flat galaxy 
density for Zspcc < 0.14. In order to evaluate the reliability of the 
identification procedure adopted in this work, we used the S3 sam- 
ple considering spectroscopic and photometric triplets detected in 
the redshift range 0.09 ^ z^pac ^0.11 (figure |5] small region) and 
0.09-CTphot < .Zphot ^ 0.11+aphot (figure|5]large region) respec- 
tively. For a discussion about the uncertainties in the photometric 
redshift see subsection 4.3. 

Using S3 sample we have obtained 64 triplet candidates with 
spectroscopic measurements. About ~ 90% of the systems had 
been recovered using photometric data. Nevertheless, there are 39 
triplets with 1 or 2 members without spectroscopic information 
due to fiber collision. Figure [8] shows two examples of photometric 
triplet candidates where one or two galaxies do not have spectro- 
scopic redshift determination. 

We have taken into account these considerations and we have 
included these objects as triplet of galaxies candidates at low red- 
shift. 



5 IMPLEMENTATION OF THE ALGORITHM AT 
INTERMEDIATE REDSHIFT (0.14 < Zphot < 0.4) 

From the discussion in Sections 4.3 and 4.4 we can see that pho- 
tometric redshifts are suitable to identify triplet of bright galax- 
ies candidates. The right panel of figure[T]show the galaxy density 
distribution as a function of redshift for the photometric S4 sam- 
ple. It can be noticed the flat behaviour of the galaxy density up to 
Zphot < 0.4. 

The unc ertainty of the phot ometric redshifts became larger at 
low redshift jO'Mill et al]|2oT l|). Taking into account these issues 
we assume a range 0.14 ^ Zphot ^ 0.4 of reliability for the deter- 
mination of photometric triplet candidates. 

We run the algorithm for S4 sample considering the restric- 
tions in AV according to the discussion of section 4.3. 

Following the steps described in section 3.2, a total of 5628 
triplet system candidates were obtained in the range 0.14 ^ 
Zphot ^ 0.4, 3771 of them with the constraints Vp < 100 h^^ kpc 
and AV < 6810km {lo-phot * c) and 1857 with the second set 
of constrains rp < 200 kpc and AV < 9310km (corre- 
sponding to ~ 1.37aphot * c). 

Figure |7] shows the distribution of the values of AVrms for 
the systems identified in the S4 sample, with the higher AV^ and 
Vp constraints (red line), as well as this distribution for systems 
identified relaxing these constraints (blue dashed line). From this 
figure it can be appreciated that both distributions show a similar 
behaviour. 



6 ISOLATION CRITERIA 

In the previous sections we have identified galaxy triplets candi- 
dates, regardless the relative spatial isolation of these systems. In 
order to build a catalogue of physical triplet systems, it is neces- 
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Figure 7. Distribution of AVrms for tlie systems identified in S4 sample. Red line: r-p < 100 h ^ kpc and AV^ < 6810km s ^ . Blue dashed line systems 
with Tp < 200 fc-i kpc and AV < 9310km s"!. 



sary to define isolation critera that ensures that the dynamics of galaxies. Nevertheless, from a dynamical point of view, there could 
these systems is not dominated by larger virialized structures where be other neighbour galaxies with lower luminosity that affects the 
these systems could be immersed. properties of the system. 



Physical triplets are a special case of a small group with three 
members. We have selected triplet candidates comprising luminous 



Hereafter, we consider two proxies for galaxy environment: 



Galaxy triplets in SDSS-DR7. 
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Figure 8. Distribution of the distance to the closest neighbour cluster for spectroscopic triplet candidates. Left panel: Distance to IZapata et al J i2009h for 
systems with 0.01 ^ Zspcc < 0.1. Right panel: Distance to GMBCG clusters for systems in the range with 0.1 ^ ^spcc < 0.14. Solid black line show the 
distribution of all the systems in that redshift range, in grey dashed line the distribution of the triplet candidates with N05 = 3 and Nl ^ 4 estimated using 
galaxies brighter than Mr = —20.5 and in solid grey line the systems with A^OS = 3 and A'^1 ^ 4 estimated using galaxies brighter than Mr = —19.5. 
Vertical black dashed line show the dc = 5 /i~^Mpc limit selected as part of the isolation criteria. 







Figure 9. Examples of tiiplets with spectroscopic and photometric measurements in the redshift range 0.01 ^ ^pho 
marked with a box and triplet galaxies members with a circle. 



^ 0.4. Galaxies with spectroscopy are 
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one associated to the local environment, and another related to 
larger scales. 

In order to describe global environment, we compute the 
distance to the closest neighbour cluster of galaxies. This quan- 
tity is related to the mass of the region where the galaxy resides 
jGonzalez & Padillal2009l) . Nevertheless, even if a triplet candidate 
system is distant from high density regions, it could still be part of a 
larger group of galaxies. To avoid such systems we have employed 
a fixed aperture method to describe the local galaxy density. 

Therefore, we have selected isolated triplets that satisfies the 
following criteria: 

(i) iVo5 = 3 

(ii) iVi s; 4 

(iii) dc ^ 5 /i"^Mpc 

and A'^i are the number of galaxies brighter than Mr = 
—20.5 within 0.5 ^~^Mpc and 1 /i^^Mpc respectively, with the 
same restrictions on /S.V used to identify triplet members, and dc 
is the distance to the closest neighbour cluster . 

These criteria were chosen from an analysis of the triplet can- 
didates at low redshifts (0.01 ^ z < 0.14) where we analysed the 
behaviour of A'os and A''i for galaxies brighter than Mr = —20.5 
and Mr — —19.5, in order to calibrate these two quantities for the 
selection of an isolation criterion to be applied at higher redshifts 
where only bright galaxies are included. In order to work with a 
complete sample of galaxies with a fainter limiting absolute mag- 
nitude than that given by the spectroscopic survey, we have con- 
sidered galaxies brighter than Mr — —19.5 up to Zspcc = 0.14, 
using both spectroscopic as well as photometric data. Accordingly, 
we estimate the A'os and TVi parameter taking into account the pho- 
tometric redshift errors. 

In order to analyse the global environment we used the pro- 
jected distance to the closest neighbour cluster in two cluster cat- 
alogues: we used a group/cluster sample extracted from the cat- 
alogue constructed by Zapata etal. (2009) updated to the SDSS- 
DR7 for 0.01 ^ ^spcc < 0.1 (dc = d:,apata)- For systems 
with 0.1 Zspoc < 0.14 we have used the Gaussian Mixture 
Brightest Cluster Galaxy (GMBCG) catalogue faao et al.l [ioTlh 
{dc — d cMBca) - 

.Zapata et al. employed a friends-of-friends algorithm 

with varying projected linking length (7, with (70 — 0.239 /i^^Mpc 
and fixed radial linking length /S.V = 450km s~^. The sam- 
ple co mprises 1514 groups and extends out to z = 0.12. Fol- 
lowing IPadilla et alj j2OI0l) . we select groups with virial masses 
Mvir > W^*h~^ Mq and define a group/c luster sample. The GM- 
BCG catalogue of clusters jHao et alj201lh comprises 55424 clus- 
ters identified using the red sequence plus Brightest Cluster Galaxy 
features, in the redshift range 0.1 ^5 2; < 0.55. The cluster abun- 
dance of GMBCG clusters is similar to IZapata et alj ^20091) pro- 
vided a mass restriction Mvir > W^'^h~^ Mq is applied to this 
sample. 

Therefore, for each triplet system candidate we compute the 
projected distance to the closest neighbour cluster considering a 
radial velocity cut AV < 1000km when spectroscopic in- 
formation is available and AV < 7000km for systems with 
photometric redshifts. 

Figure [8] shows the distribution of the distance to the closest 
neighbour cluster (left Zapata et al. 2009 and right Hao et al. 20 1 1 ), 
for the spectroscopic triplet candidates. In black solid line the dis- 
tribution for systems in the redshift range 0.01 ^ 2:spec < 0.1 (left 
panel) and 0.1 ^ ^spcc < 0.14 (right panel). In these figures, grey 
dashed lines represent the distribution of triplet candidates that sat- 



isfies A'os ~ 3 and A^i ^ 4, estimated using galaxies brighter 
than Mr — —20.5. The distribution for systems with = 3 
and A'^i ^ 4, estimated with galaxies A/,. ^ —19.5 is represented 
in grey solid line. In both panels of this figure it can be appreci- 
ated that the limit dc — 5 h^^Mpc (vertical black dashed line) 
encloses most of the systems with A'os ~ 3 and A''i ^ 4 (galaxies 
Mr ^ -19.5 as well as galaxies with Mr ^ -20.5). 

We have considered this limit in order to select a safe isola- 
tion criterion, for triplet candidates at higher redshift (0.14 5^ z ^ 
0.4), where, by completeness, the selection of isolated triplets with 
A'os = 3 and A'^i ^ 4 can be estimated only by using galaxies 
brighter than Mr = -20.5. 

Taking a distance grater than 5 /i~^Mpc, we are ensuring that 
our triplets do not belong to any system of groups or clusters. 
On the other hand, by restricting the number of galaxies within 
0.5 /i~^Mpc and 1 /i~^Mpc leave us with a system composed by 
only three luminous galaxies. 

By imposing the isolation criterion to the sample of triplet of 
galaxies candidates we obtain 1092 triplets in our final catalogue, 
95 triplets in the redshift range 0.01 ^ Zspcc < 0.14 and 997 in 
the range 0.14 ^ 2:phot 0.4. Figure |9] shows three examples of 
the systems at low and intermediate redshifts. From these figures 
it can be appreciated the high degree of isolation of the systems. 
In a forthcoming paper (paper II, Duplancic et. al in preparation), 
we will analyse different spectro-photometric properties of these 
isolated triplet of galaxies systems. 

In the Appendix we summarise some properties of 20 triplets 
from the final catalogue derived in this work. Columns 1 and 2 right 
ascension and declination (J2000) galaxy positions in hexadecimal 
format; column 3 redshift of the central galaxy, and column 4 num- 
ber of triplet galaxy members with spectroscopic measurements 
The full catalog is available electronically. 



7 RESULTS AND DISCUSSION 

In this work we have developed algorithms for the identification 
of galaxy triplet candidates in the Sloan Digital Sky Survey Data 
Release Seven. 

Taking into account the completeness of the spectroscopic 
data in the range 0.01 ^ ^spcc < 0.14 down to absolute magni- 
tude Adr = —20.5, we identified 273 triplet candidates using only 
the spectroscopic data. 

We have considered the incompleteness of these spectroscopic 
triplets due to fiber collision by using photometric as well as spec- 
troscopic data. We find that 90% of the spectroscopic systems can 
be recovered using the photometric data and, additionally, 39 new 
triplet candidates with 1 or 2 members without spectroscopic infor- 
mation were included into this sample. 

In order to explore the completeness and contamination of our 
procedure, we have employed a mock catalogue with similar vol- 
ume than the spectroscopic SI sample. We find a high level of com- 
pleteness (~ 80%) and low contamination (~ 5%) which gives 
confidence to the methods adopted and the results obtained. 

Using photometric information we have extended the sample 
of triplet candidates in the spectroscopic redshift range up to sig- 
nificantly larger depths (0.14 ^ Zphot ^ 0.4). Within this larger 
redshift range we identified 5628 triple system candidates. 

In order to build a catalogue of triplets of galaxies which fulfil 
the main features of physical systems, we defined an isolation cri- 
terion based on both, the distance to the closest neighbour cluster 
and a fixed aperture local density estimate. This criterion ensures 



that the dynamics of the systems is not dominated by larger virial- 
ized structures and also avoids triplets contaminated by bright close 
neighbour galaxies. 

By imposing the isolation criterion to the sample of triplet of 
galaxies candidates we obtain 1092 triplets in our final catalogue, 
95 triplets in the redshift range 0.01 ^ Zspec < 0.14 and 997 in 
the range 0.14 ^ Zphot ^ 0.4. The results can be appreciated by 
inspection of Figure (TO) which shows the spatial number density 
of isolated triplets as a function of redshift, for the low redshift 
sample. Red dots corresponds to the spectroscopic triplets and blue 
dots to the spatial number density with the new triplets included 
after the fiber collisions analysis (see section 4.4). From this figure 
it can be seen that photometric redshift information is useful in 
order to complete the sample of triplets at low redshift. 

Figurefm shows the spatial number density of isolated triplets 
as a function of redshift for the intermediate redshift sample. By 
comparison to figure[TO]it can be appreciated that the spatial num- 
ber density of both, the intermediate and the low redshift photomet- 
ric triplet samples, smoothly match at z ^ 0.14. The observed flat 
behavior of the number density up to jZphot ~ 0.4 shows that our 
methods are adequate to identify triplets reliably without important 
systematic effects. 

Our results show that photometric redshifts provide very use- 
ful information, allowing to asses fiber collision incompleteness 
and extend the detection of triplets to large distances where spec- 
troscopic redshifts are not available. 

The final catalogue comprise 1092 triplets of galaxies in the 
redshift range 0.01 ^ z 0.4. This sample is suitable for statisti- 
cal analysis, which will be performed in a forthcoming paper, and 
may serve as targets for follow up observations. 
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Figure 10. Spatial number density of isolated triplet of galaxy at low redshift. Red dots represents systems with spectroscopic information and in blue the 
spatial number density with the systems incorporated by using photometric information. The eiTor bars represent Icr uncertainties calculated using bootstrap 
techniques. Vertical lines shows the redshift range analysed in section 4.4. 




Figure 11. Spatial number density of isolated triplet of galaxy at intermediate redshift. The error bars represent Icr uncertainties calculated using bootstrap 
techniques. 
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Table Al. Examples of triplets of galaxies. Column 1 and 2: right ascension 
and declination galaxy position. Column 3: redshift of the central galaxy. 
Column 4: number of triplet galaxy members with spectroscopic measure- 
ments. 
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APPENDIX A: CATALOGUE SUMMARISE 



We present a short table (table IaH with some properties of 20 
triplets from the final catalogue derived in this work. 



